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Abstract 

The final results of the MACRO experiment on atmospheric neutrino oscillations are 
\ presented and discussed. The data concern different event topologies with average neutrino 

£NJ ■ energies of ~ 3 and ~ 50 GcV. Multiple Coulomb Scattering of the high energy muons in 

absorbers was used to estimate the neutrino energy of each event. The angular distributions, 
the L/E u distribution, the particle ratios and the absolute fluxes all favour — > v T oscil- 
lations with maximal mixing and Am 2 = 0.0023 eV 2 . A discussion is made on the Monte 
■ Carlos used for the atmospheric neutrino flux. Some results on neutrino astrophysics are 

""^1 \ also briefly discussed. 

Sjj ■ Keywords: Astroparticle physics; underground detectors; neutrino oscillations. 
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D 

1 Introduction 

> ■ 

^ . MACRO was a large area multipurpose underground detector designed to search for rare 

events and rare phenomena in the penetrating cosmic radiation. The experiment obtained 
important results on atmospheric oscillations and performed £/„ astronomy studies, indirect 
searches for WIMPs and searches for low energy v e from stellar gravitational collapses 0. 

The detector was located in Hall B of the underground Gran Sasso Lab at an average rock 
overburden of 3700 m.w.e.; it started data taking with part of the apparatus in 1989; it was 
completed in 1995 and was running in its final configuration until December 2000. 

The detector had global dimensions of 76.6 x 12 x 9.3 m 3 and provided a total acceptance to 
an isotropic flux of particles of ~ 10, 000 m 2 sr. The detector was composed of three horizontal 
layers of liquid scintillation counters, 14 layers of limited streamer tubes and one layer of nuclear 
track detectors, Fig. ^ Vertically it was divided into two parts: the lower part contained 10 
layers of streamer tubes, 7 layers of rock absorbers and 2 layers of liquid scintillators; the upper 
part was empty, contained the electronics and was covered by 1 layer of scintillators and 4 layers 
of streamer tubes. The sides of the detector were covered with 1 layer of scintillators and 6 layers 
of limited streamer tubes, so as to obtain a closed box structure. Each of the subdetectors could 
be used in "stand-alone" and in "combined" mode. 



*see Ref. [1] for a list of MACRO Authors and Institutions 
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In the following we shall briefly recall neutrino oscillation formulae and atmospheric neutri- 
nos; then we shall discuss our main results including the ratio of vertical to horizontal muons 
and the study of the L/E v distribution. After a detailed analysis of the present situation of the 
different Monte Carlos (MCs), we summarize our atmospheric neutrino oscillation results and 
we conclude with a short discussion of some neutrino astrophysics studies. 




Figure 1: Cross section of the detector, which also shows the various atmospheric neutrino topologies 
measured. The horizontal scintillators are indicated as dashed boxes, the streamer tubes as horizontal 
thick black lines, the rock absorbers as boxes with horizontal thin lines. The sides of the detector were 
formed by 1 layer of scintillators and 6 layers of streamer tubes. 



2 Atmospheric neutrino oscillations. Monte Carlos 

Neutrino oscillations. If neutrinos have non-zero masses, one has to consider the 3 weak 
flavour eigenstates u e , v^, u T and the 3 mass eigenstates u\, V2, ^3- The flavour eigenstates 
are linear combinations of the mass eigenstates v m . Neutrino oscillations depend on six 
independent parameters: two mass-squared differences, Amj 2 and Am^, three mixing angles 
#12, #13, #23 and the CP-violating phase 5. In the simple case of two flavour eigenstates (z^, v T ) 
which oscillate with two mass eigenstates (f2, ^3), and 5 = one has: 

Vjx = V2 COS 6>23 + Vj, Sin 6*23 (1) 
v T = -u 2 sin #23 + ^3 cos 6>23 (2) 
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The survival probability of a v„ beam is 

P(y» -+Vy)~\- sin 2 29 23 sin 2 I ^ J (3) 

where Am^ = m§ — to 2 , and L is the distance travelled by the neutrino of energy E u from 
production to detection. The two-neutrino approximation is adequate for discussing our data. 

Atmospheric neutrinos. High energy (HE) primary Cosmic Rays (CRs), protons and 
nuclei, interact in the upper atmosphere producing a large number of pions and kaons, which 
decay yielding muons and muon neutrinos; then muons decay yielding Un's and zVs. The ratios 
of their numbers are N Vii /N Ve ~ 2 and N u /Nj7 ~ 1. The atmospheric neutrinos are produced in 
a spherical shell at about 10-20 km above ground and proceed towards the Earth. Atmospheric 
neutrinos have energies from a fraction of GeV up to more than 100 GeV and they travel distances 
L from few tens of km up to ~ 13000 km. An underground detector is "illuminated" by a flux 
of neutrinos from all directions and it can make oscillation studies for 1 < L/E v < 10 4 km/GeV. 

MACRO detected upgoing muon neutrinos via charged current interactions, — » /i; the 
upgoing muons were identified with the streamer tube system (for tracking) and the scintilla- 
tor system (for time-of-flight measurement). The events measured and expected for the three 
measured topologies, indicated in Fig. ^and the L/E v distribution deviate from Monte Carlo 
expectations without oscillations, Figs El 0] and Tables 1, 2; the deviations point to the same 
Vfi — ► v T oscillation scenario [2]-[H]- 

Monte Carlos. The measured data of Fig. EJwere compared with different MC simulations. 
In the past we used the neutrino flux computed by the Bartol96 group [0] and the GRV94 ^0] 
parton distribution. For the low energy channels the cross sections in were used; the 
propagation of muons to the detector was done using the energy loss calculation by Lohmann 
et al. ^2]- The total systematic uncertainty in the predicted flux of upthroughgoing muons, 
adding in quadrature the errors, was estimated to be ±17 %; this is mainly a scale error that 
does not change the shape of the angular distribution. The error on the shape of the distribution 
is ~ 6%. The detector was simulated using GEANT ^H]- A similar MC (Honda96) was used by 
the Superkamiokande Collaboration |14j . 

Recently new improved MC predictions were made available by the Honda Jl] and Fluka 
|15j groups. They include three dimensional calculations of hadron production and decays 
and of neutrino interactions, improved hadronic model and new fits of the primary cosmic ray 
flux. The two MCs yield predictions for the non oscillated and oscillated Un fluxes in complete 
agreement with each other, equal to within few %, see Figs. |2h> b. The shapes of the angular 
distributions for oscillated and non oscillated Bartol96, new Fluka and new Honda fluxes are the 
same to within few %, Fig. b. The absolute values of our upthroughgoing muon data are 
about 20 — 30% above those predicted by the new Fluka and Honda MCs; this situation is also 
true for the new Superkamiokande data |16j : the measurements of the flux value at cosi? = 
performed by various experiments do not seem to be conclusive |17j . Note in Fig. |2K that the 
new Fluka muon flux with an older cosmic ray fit is considerably above the new Fluka flux with 
the new CR fit. Our high energy data suggest that the new Honda and Fluka predictions 
should be raised by about 20 — 25%, probably because of the used new CR fit. The inclusion of 
the new ATTIC Collaboration measurements of primary CRs seems to lead to the old energy 
dependence of E~ 2 - 71 [TBI - Thus the Bartol96 MC may probably still be used for the prediction of 
the absolute flux, besides predicting well the shape of the angular distribution. In the following 
we shall use mainly the predictions of this MC. It should be noted that the evidence for neutrino 
oscillations rests primarely with the shape of the angular distribution and this is the same in all 
MC calculations. 
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3 MACRO results on atmospheric neutrinos. 



The upthroughgoing muons come from interactions in the rock below the detector; the Vu's 
have a median energy E v ~ 50 GeV; muons with > 1 GeV cross the whole detector. The 
data taking livetime was slightly over 6 years (full detector equivalent). The data, Table 1 and 
Fig. 12^,, deviate in shape and in absolute value from the Bartol96 MC predictions (see Section 
2). This was first pointed by MACRO in 1995 gj. 
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Figure 2: Zenith distributions for the MACRO data (black points) for (a) upthroughgoing, (b) semicon- 
tained and (c) up-stopping muons + down semicontained. The dashed lines are the no-oscillation MC 
predictions Bartol96 in (a) and new FLUKA in (b) and (c) (with scale error bands); the solid lines refer 
to — » v T oscillations with maximal mixing and Am 2 = 2.3 • 10 -3 eV 2 . 

A large number of possible systematic effects and backgrounds that could affect the measure- 
ments were studied. One of the most significant checks was performed using only the scintillator 
system with two different and independent electronic systems (HI IE]- Two different final analyses 
of the upthrough data have been performed: they agree to within 5%. 

Assuming no oscillations, the number of expected upthroughgoing muon events integrated 
over all zenith angles from Bartol96 without oscillations is 1169; the measured number is 857, 
Table 1. Thus the ratio of the observed number of events to the Bartol96 expectation is 0.73. 

f„ — > v T versus i/„ — > v s . Matter effects due to the difference between the weak interaction 
effective potential for muon neutrinos with respect to sterile neutrinos (u s ) would produce a 
different total number and a different zenith angle distribution of upthroughgoing muons. The 
ratio Ri = Vertical/Horizontal = 2V(— 1 < cosO < -0.7)/iV(-0.4 < cosQ < 0) was used to 
test the f/j, — > v s oscillation hypothesis versus — ► v T [2], [S], |Hj- The measured value is 
Rmeas = 1-38; it should be compared to i?™ n = 1.61 and R™™ii e = 2.03, which are the minimum 
values of the ratios for Un — ► v T an d Vn — > ^sterile oscillations for Am^ = 2.3 • 10~ 3 eV 2 and 
maximal mixing. For — > v T , the one-sided probability to measure a value lower than R mea s 
is 7.2%. For — > v s the probability is ~ 1.5 • 10~ 4 ; P^ est is about 480 times larger than 
Pster- Therefore i/„ — > v s oscillations (with any mixing) are excluded at about 99.8% c.l. with 
respect to Vn — > v T oscillations with maximal mixing 8 u . 

Oscillation probability as a function of the ratio L/E u . E u was estimated by mea- 
suring the muon energy, En, by means of the muon Multiple Coulomb Scattering (MCS) in 
the rock absorbers in the lower MACRO. The best method employed the streamer tubes in 
"drift mode", using the special electronics designed to search for magnetic monopoles ^Hj. The 
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(a) (b) 

Figure 3: (a) Comparison of the measured angular distribution for upthroughgoing muons and the 
oscillated flux obtained from different MCs (for Am 2 — 2.3 • 10 -3 eV 2 and maximal mixing). The solid 
line is the Bartol96 flux, the dotted and the dashed lines at the bottom are the new Fluka and Honda 
fluxes; the dotted-dashed line in the middle is the Fluka flux using a different (older) fit to the cosmic 
ray data, (b) Comparison of our measurements with the Bartol96 and the new Honda 2001 oscillated 
and non oscillated fluxes. 



Topology 


Measured events 


{Data /MC) noosc 


Upthroughgoing muons 


857 


0.73 


Semicontained (UP) 


157 


0.55 


Upstopping + Down 


262 


0.74 


semicontained (ID + UGS) 







Table 1: For each topology the table gives the number of events measured and the ratios (Data/M C) noosc 
using the Bartol96 MC. 

space resolution achieved was ~ 3 mm. For each muon, seven variables were given in input 
to a Neural Network (NN) previously trained to estimate muon energies with MC events of 
known input energy crossing the detector at different zenith angles. The distribution of the 
ratio R = (Data/MC noosc ) obtained by this analysis is plotted in Fig. 0]as a function of (L/E v ) 
7 . Notice that the data extend from (L/E„) ~ 30 km/GeV to 5000 km/GeV. 

The Internal Upgoing (IU) muons (~ 157 events, 285 expected by Bartol96 and 235 expected 
by new FLUKA, Table 1) come from ~ 4 GeV u^s interacting in the lower apparatus, Fig. ^ 
Compared to the no-oscillation prediction there is a reduction of about a factor of two in the 
flux of these events, without any appreciable distortion in the shape of the zenith distribution, 
Fig. EJd 5 , 8 . The MC predictions for no oscillation in Figs. 03 and[2j: are given by the dashed 
lines with a 21% systematic scale band. Notice that at these energies, the new Honda and Fluka 
MC predictions are low by about 12%. 

The upstopping muons (UGS) are due to ~ 3 GeV z^'s interacting below the detector and 
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Figure 4: Ratio (Data/MC Bartol96) as a function of estimated L/E v for the upthrougoing muon sample 
(black circles) and the semicontained up-/i (open circle). For upthroughgoing muons the muon energy 
was estimated by multiple Coulomb scattering and E v by MC methods. The shaded regions represent 
the uncertainties in the MC predictions assuming sin 2 29 = 1 and Am 2 — 0.0023 eV 2 . The horizontal 
dashed line at Data/MC=l is the expectation for no oscillations. 



yielding upgoing muons stopping in the detector. The semicontained downgoing muons (ID) are 
due to z^-induced downgoing muon tracks with vertex in the lower MACRO. The two types of 
events (262 compared to 354 expected by Bartol96 and 315 by FLUKA, Table 1) are identified 
by means of topological criteria; the lack of time information prevents to distinguish the two 
sub-samples. The upgoing v^'s should all have oscillated completely, while the downgoing z/« do 
not. The zenith distribution shows, as expected, a uniform deficit of about 25 % of the measured 
number of events with respect to the no-oscillation prediction, Fig. [2k 0, jSj. 

4 Determination of the oscillation parameters 

In the past, in order to determine the oscillation parameters, MACRO made fits to the shape of 
the upthroughgoing muon distribution and to their absolute flux compared to the Bartol96 MC 



Ratio 


Rmeas ± 0~stat 


i?0 


< Rmeas probab 
a priori II 


Vertical /Horizontal 


1.38 ±0.12 


2.11 


6.4 x 10~ e 


Nloui/Nfiigh 


0.85 ±0.16 


1.5 


7.7 x 10~ 3 


InUp/ '(InDown + UpStop) 


0.60 ± 0.06 


0.745 


3.1 x 10~ 2 


Combination 






3.4 x 10~ 7 



Table 2: Ratios between different event categories: i? meas is the measured value, i?o is the expected 
value from the Bartol96 MC without oscillations; using the new FLUKA MC the ratios differ by at most 
5%. The last column gives the one-sided probability compatible with a statistical fluctuation computed 
according to the a priori assumptions explained in the text. 
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Figure 5: Interpolated qualitative 90% C.L. contour plots of the allowed regions in the Am 2 — sin 2 29 
plane (a) for the MACRO data using only the ratios R\,R%, R3 (outer continuous line) and using also the 
absolute values assuming the validity of the Bartol96 fluxes, (b) Comparison of the allowed regions for 
SK, MACRO and Soudan 2 (All contours have been smoothed out and must be considered as qualitative). 



prediction (with a 17% systematic scale uncertainty). The other data were only used to verify 
the consistency and make checks. The result was Am 2 = 0.0025 eV 2 and maximal mixing. 

Later, in Ej the high energy data for the zenith distribution, Ri = (N vert / Nh or iz) , were 
combined with the HE neutrino energy measurement, R 2 = Ni^/N^igh- They are independent 
of the neutrino absolute flux; with these two informations the significance of observation of 
neutrino oscillations is about 4.7 a. 

In order to reduce the effects of possible uncertainties in the MCs (to about 6%) we now use 
the following three independent ratios |Hj. 

(i) High Energy Data: zenith distribution ratio: R\ = N ver t/Ni lor 

(ii) High Energy Data, neutrino energy measurement ratio: R2 = Ni ow /Nhi g h 
(hi) Low Energy Data: Ratio R3 = (Data/MC)iu/(Data/MC)iD+UGS- 

Table 2 gives for each ratio the measured value, R meas , the prediction, Rq, for no os- 
cillations using either the Bartol96 or the new FLUKA MC; the probability value for no 
oscillations obtained allowing the predicted ratio to fluctuate around the mean value (con- 
sidering the total number of events equal to the predicted value without oscillations). The 
no oscillations hypothesis has a probability P ~ 3.4 x 10~ 7 (Table 2), ruling out the no- 
oscillation hypothesis by ~ 5a. The formula used for combining independent probabilities is 
P = PiP 2 P 3 (l - In Pi P2P3 + l/2(lnP 1 P 2 P 3 ) 2 ) m- 

By fitting the three ratios to the — > v T oscillation formulae we obtain sin 2 2$ = 1, Am 2 = 
2.3 • 10 -3 eV 2 and the allowed region indicated in Fig. [5^. 

We could also add the information on the absolute flux values of the 

(iv) high energy data, assuming the validity of the Bartol96 flux and a systematic scale error 
of 17%, R 4 = N meas /N Bartol _ nooscill ~ 0.73, P 4 ~ 0.07 

(v) low energy semicontained muons, with a systematic scale error of 21%, 

R5 = N meas /N Bartol _ noosciU ~ 0.67, P 5 ~ 0.08 (At these low energies the Bartol96, Honda 
2001 and Fluka 2002 MC neutrino fluxes are essentially equal). 
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Fig. shows the MACRO allowed regions for the v„ — > v T oscillation parameters in the 
Am 2 — sin 2 26 plane using only the three ratios R\, R2 and i?3 (continuous line) and using also 
the two estimated absolute values, assuming the validity of the Bartol96 fluxes (dashed line). 
In Fig. \5}p are compared the MACRO, SuperKamiokande (SK) ^f)j and Soudan 2 (2^ allowed 
regions. All limiting lines are qualitative smoothed interpolations. 

5 Neutrino Astrophysics with MACRO 

Search for Astrophysical Sources of High Energy Muon Neutrinos. High energy v^s 
are expected to come from several galactic and extragalactic sources. Neutrino production re- 
quires astrophysical accelerators of protons and astrophysical "beam dumps". A sensitive search 
was made for upgoing muons produced by neutrinos coming from celestial sources, interacting 
below the detector. 90% c.l. upper limits were established on the muon fluxes from specific 
celestial sources; the limits are in the range cm s . In the case of GX339-4 

(a = 255.71°, 5 = -48.79°) and Cir X-l (a = 230.17°, 5 = -57.17°), we have 7 events: they 
have been considered as background, therefore the upper flux limits are higher; but the events 
could also be indication of signals |221 12*1 13]. The pointing capability of the detector was tested 
by determining the shadows of the Moon and of the Sun on the primary cosmic ray flux |23| 

A search for time coincidences of the upgoing muons with 7-ray bursts was made. No statis- 
tically significant time correlation was found. A search was also made for a diffuse astrophysical 
neutrino flux for which an upper limit at the level of 1.5 • 10~ 14 cm _2 s _1 was established |24| . 

Indirect Searches for WIMPs. Weakly Interacting Massive Particles (WIMPs) could be 
part of the galactic dark matter; they could be intercepted by celestial bodies, slowed down and 
trapped in their centers, where WIMPs and anti- WIMPs could annihilate and yield neutrinos 
of multi GeV energy, in small angular windows from their centers. One of the best WIMP 
candidate could be the lowest mass neutralino. 

A search was made for Vu — ► /i coming from the center of the Earth using 10° — 15° cones 
around the vertical; we obtained muon flux limits of ~ 10 -14 cm -2 s _1 . The limits were compared 
with the predictions of a supersymmetric model by Bottino et al.; our data eliminate a sizable 
range of parameters used in the model [23 Ej • A similar procedure was used to search for muon 
neutrinos from the Sun: the upper limits are at the level of about 1.5 • 10~ 14 cm -2 s^ 1 |25| . 

Neutrinos from Stellar Gravitational Collapses. A stellar gravitational collapse of the 
core of a massive star is expected to yield a burst of all types of neutrinos and antineutrinos with 
energies of 5 — 30 MeV in ~ 10 s. The v e 's can be detected via v e + p — > n + e + in the liquid 
scintillator. About 120 v e events could have been detected in our 580 t liquid scintillator for a 
stellar collapse at the center of our Galaxy. Two separate electronic systems were used. Both 
had an energy threshold of ~ 7 MeV and recorded pulse shape, charge and timing informations 
of the positron. Following a > 7 MeV trigger, the first system PHRASE lowered its threshold to 
~ 1 MeV for 800 [is in order to detect (with a ~ 25% efficiency) the 2.2 MeV 7 released in the 
capture reaction n + p — > d + 7 induced by the neutron produced in the primary process. A 
redundant supernova alarm system was in operation, alerting the physicists on shift; a procedure 
was defined to alert the various supernova observatories around the world MACRO was 
completed at the end of 1994; the first parts of the detector started operation in 1989. No stellar 
gravitational collapses in our Galaxy were observed from 1989 to 2000. 
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6 Conclusions 



The MACRO detector took data from 1989 till the end of year 2000. Important results were 
obtained in all the items listed in the proposal. 

In particular for atmospheric neutrino oscillations we analyzed different event topologies, 
different energies, exploited muon Coulomb multiple scattering in the detector to measure muon 
energies and studied the effects of using different MCs. All data are in agreement with the 
hypothesis of — » v T oscillations, with maximal mixing and Am^ ~ 2.3 • 10~ 3 eV 2 , and rule 
out the no oscillation hypothesis by ~ 5a. 

Studies were made on muon neutrino astronomy |22j . |23j . seasonal, solar and sidereal vari- 
ations [5] and on possible Dark Matter candidates : WIMPs (looking for muon neutri- 
nos from the centers of the Earth and of the Sun), Nuclearites and Q-balls [2J, |27j. In the 
search for GUT Magnetic Monopoles MACRO obtained the best existing direct flux upper limit 
($ < 1.4 x 10 -16 cm -2 s -1 sr -1 ) over the widest range (4 x 10~ 5 < /3 < 1) [T^] . 
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